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Given a positive integer n, let P(n) denote the largest prime factor of n and S/n) denote the 
smallest integer m such that n divides m! 


This paper extends earlier work [1] on the average value of the Smarandache function S/n) 
and is based on a recent asymptotic result [2]: 
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due to Ford. We first prove: 
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where ¢ (x) is the Riemann zeta function. In particular, 


. In(N) m? 
in = “= 082246703... 
lim = - E(S(N)) = = = 082246703 
In(V 3 
lim D, Var(S(N)) = £2) = 0.40068563... 
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Sketch of Proof. On one hand, 
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The above summation, on the other hand, breaks into two parts: 
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The second part vanishes: 
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P(n)<S(n) 








while the first part is bounded from above: 
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A formula for L(k) was found by Knuth and Trabb Pardo [3] and the remaining second- 
order details follow similarly. 


Observe that the ratio /Var(S(N)) / E(S(N)) > 2% as N >, which indicates that the 
traditional sample moments are unsuitable for estimating the probability distribution of 
S(N). An alternative estimate involves the relative number of digits in the output of S per 
digit in the input. A proof of the following is similar to [1]; the integral formulas were 
discovered by Shepp and Lloyd [4]. 


Theorem 2. 
0.62432998 if 


k=1 
nS) (xk eho 0.42669576 if k=2 

lim n(S(V)) |_ 2 hl oe 2p ees esiseiel- i as 
i : 0.24387660 if k=4 

0.19792289 if k=5 


The mean output of S hence has asymptotically 62.43% of the number of digits of the 
input, with a standard deviation of 19.21%. A web-based essay on the Golomb-Dickman 
constant 0.62432998... appears in [5] and has further extensions and references. 
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